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Abstract 

This paper describes a study of B meson decays to the pseudoscalar-vector final 
state pir using 31.9 x 10 6 BB events collected with the Belle detector at KEKB. The 
branching fractions B(B + -► p°ir + ) = (8.0±|;g±g;^) x 1(T 6 and B(B° -» p ± ^) = 
(20.8 

ies!3^3!i) x 10 6 are obtained. In addition, a 90% confidence level upper limit 
of B(B° -► p°7r°) < 5.3 x lO" 6 is reported. 

Key words: pn, branching fraction 
PACS: 13.25.hw, 14.40.Nd 



The decays of _B-mesons to pseudoscalar and vector particles provide op- 
portunities for investigating the phenomenon of CP violation. Of particular 
interest are the quasi-two-body B — > pit decays to three-pion final states. 
These can be used to measure the angles 02 and 03 of the unitarity tri- 
angle: 02 can be determined from a full Dalitz plot analysis of the modes 
B° — > p ±r K T and B° — > p°7r° [1]; 03 can be extracted from the interference 
of B + — > Xco7T + and B + — > p°n + in B + — > 7r + 7r _ 7r + decays [2]. While the 
data set used here is too small for this type of analysis, it can be used to 
study the branching fractions of these decays. The ratio of branching fractions 
R = B(B° — > p ± 7i T )/B(B + — > p°7r + ) is particularly interesting as theoreti- 
cal predictions for the value of R vary over a wide range depending on the 
assumptions made in its calculation [3-8]. 

In this paper, a search for B meson decays of the type B — > pit is described. 
Both the neutral (B° — > p ±r n T and B° — > p°7r°) and charged (£? + — > p°7r + ) 
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modes are examined. Here and throughout the text, inclusion of charge con- 
jugate modes is implied and for the neutral decay, B° — > p ± 7r T , the notation 
implies a sum over both the modes. The data sample used in this analysis 
was taken by the Belle detector [9] at KEKB [10], an asymmetric storage 
ring that collides 8 GeV electrons against 3.5 GeV positrons. This produces 
T(4iS) mesons that decay into B°B° or B + B~ pairs. The data sample has an 
integrated luminosity of 29.4 fb _1 and consists of 31.9 x 10 6 BB pairs. 

The Belle detector is a general purpose spectrometer based on a 1.5 T su- 
perconducting solenoid magnet. Charged particle tracking is achieved with a 
three-layer double-sided silicon vertex detector (SVD) surrounded by a central 
drift chamber (CDC) that consists of 50 layers segmented into 6 axial and 5 
stereo super-layers. The CDC covers the polar angle range between 17° and 
150° in the laboratory frame, which corresponds to 92% of the full centre of 
mass (CM) frame solid angle. Together with the SVD, a transverse momen- 
tum resolution of (<J pt /pt) 2 = (0.0019pt) 2 + (0.0030) 2 is achieved, where p t is 
in GeV/c. 

Charged hadron identification is provided by a combination of three devices: 
a system of 1188 aerogel Cerenkov counters (ACC) covering the momentum 
range 1-3.5 GeV/c, a time-of-flight scintillation counter system (TOF) for 
track momenta below 1.5 GeV/c, and dE/dx information from the CDC for 
particles with very low or high momenta. Information from these three devices 
is combined to give the likelihood of a particle being a kaon, L^, or pion, 
L n . Kaon-pion separation is then accomplished based on the likelihood ratio 
L n / (L^+Lk)- Particles with a likelihood ratio greater than 0.6 are identified as 
pions. The pion identification efficiencies are measured using a high momentum 
D* + data sample, where D* + — > D°tt + and D° — > K~tt + . With this pion 
selection criterion, the typical efficiency for identifying pions in the momentum 
region 0.5 GeV/c < p < 4 GeV/c is (88.5 ±0.1)%. By comparing the D* + data 
sample with a Monte Carlo (MC) sample, the systematic error in the particle 
identification (PID) is estimated to be 1.4% for the mode with three charged 
tracks and 0.9% for the modes with two. 

Surrounding the charged PID devices is an electromagnetic calorimeter (ECL) 
consisting of 8736 CsI(Tl) crystals with a typical cross-section of 5.5 x 5.5 cm 2 
at the front surface and 16.2 Xo in depth. The ECL provides a photon energy 
resolution of (a E /E) 2 = 0.013 2 + (0.0007/E) 2 + (0.008/E 1 / 4 ) 2 , where E is in 
GeV. 

Electron identification is achieved by using a combination of dE/dx measure- 
ments in the CDC, the response of the ACC and the position and shape of 
the electromagnetic shower from the ECL. Further information is obtained 
from the ratio of the total energy registered in the calorimeter to the particle 
momentum, E/pi ab . 
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Charged tracks are required to come from the interaction point and have 
transverse momenta above 100 MeV/c. Tracks consistent with being an elec- 
tron are rejected and the remaining tracks must satisfy the pion identification 
requirement. The performance of the charged track reconstruction is studied 
using high momentum r] — > 77 and r] — > 7r + 7r~7r° decays. Based on the relative 
yields between data and MC, we assign a systematic error of 2% to the single 
track reconstruction efficiency 

Neutral pion candidates are detected with the ECL via their decay it — > 77. 
The 7T° mass resolution, which is asymmetric and varies slowly with the n° 
energy, averages to a = 4.9 MeV/c 2 . The neutral pion candidates are selected 
from 77 pairs by requiring that their invariant mass to be within 3a of the 
nominal ir° mass. 

To reduce combinatorial background, a selection criteria is applied to the pho- 
ton energies and the 7r° momenta. Photons in the barrel region are required to 
have energies over 50 MeV, while a 100 MeV requirement is made for photons 
in the end-cap region. The n° candidates are required to have a momentum 
greater than 200 MeV/c in the laboratory frame. For 7T°s from BB events, the 
efficiency of these requirements is about 40%. Furthermore, the 7r° helicity, de- 
fined as the absolute cosine of the angle between the photon direction in the 
7T° rest frame and the 7r° direction in the lab frame, is required to be less than 
0.95. The behavior of the n reconstruction efficiency as a function of momen- 
tum is studied using D* + — > D°ir + decays by comparing the relative yields 
between data and MC in the D° — > K~it + it and K~ir + sub-channels. The 
residual systematic uncertainties are 9.8% for low momentum 7r°s (p < 500 
MeV/c) and 7.7% otherwise. 

Neutral p candidates are reconstructed using the decay p° — > it + it~ with the 
invariant mass of the charged pion pair required to be between 0.6 GeV/c 2 
and 0.95 GeV/c 2 . The decay p + — > tt°-k + is used to reconstruct charged p 
candidates with the invariant mass of the pion pair required to be between 
0.62 GeV/c 2 and 0.92 GeV/c 2 . 

Candidate B mesons are identified by pairing p and n candidates that pass the 
above criteria. The candidates are then selected using the beam-constrained 
mass Mb c = \J E^ clim — p 2 B and the energy difference AE = E B — E hciim . Here, 
Pb and E B are the momentum and energy of a B candidate in the CM frame 
and -Ebeam is the CM beam energy. An incorrect mass hypothesis for a pion or 
kaon produces a shift of about 46 MeV in AE, providing extra discrimination 
between these particles. The width of the M bc distributions is primarily due 
to the beam energy spread and is well modelled with a Gaussian of width 
3.3 MeV/c 2 for the modes with a neutral pion and 2.7 MeV/c 2 for the mode 
without. The AE distribution is found to be asymmetric with a small tail 
on the lower side for the modes with a 7r°. This is due to 7 interactions with 
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material in front of the calorimeter and shower leakage out of the calorimeter. 
The AE distribution can be well modelled with a Gaussian when no neutral 
particles are present. Events with 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 and \AE\ < 
0.3 GeV are selected for the final analysis. 

The dominant background comes from continuum e + e~ — > qq (q — u, d, s, 
c) production. This background is suppressed using variables that quantify 
the difference in the event topology between spherical BB events and jet-like 
qq events. The most powerful suppression is achieved using variables derived 
from the Fox- Wolfram moments [11] but altered to contain the additional 
information of whether or not a track comes from the B candidate in the 
event [12]. These modified moments enhance the discriminating power of the 
original Fox- Wolfram moments and are defined as 

R = J2i,k \pi\\pk\Pi(co8 9ik) = Eij Hbj|P;( cos %) 
J2i,k\Pi\\Pk\ ' 1 Ei,j\pi\\pj\ 

where p indicates the particle momentum, / runs from 1 to 4 and Pi is the Leg- 
endre polynomial of /th order, k runs over the daughters of a B candidate and 
% and j enumerate photons and charged particles in the event not associated 
with the B candidate. Since Ri, R 3 and r\ are found to be correlated with 
Mb c , they are not used. Another event topology variable used is S±, defined as 
the scalar sum of the transverse momenta (with respect to the thrust axis) of 
the particles outside a 45° cone around the thrust axis, divided by the scalar 
sum of the momenta of all the particles. These variables then are combined to 
form a Fisher discriminant [13] 

T= E a i R i+ E Piri + lS ± , 

1=2,4 1=2,3,4 

where ai, f3i and 7 are optimized to maximize the separation between signal 
and continuum events. Also used for continuum suppression is the B flight 
direction, cos 9b, defined as the cosine of the angle the B meson momentum 
vector makes with the positron beam axis. As cos 9b is uncorrelated with the 
other continuum suppression variables, it is left out of T . 

Probability density functions (PDFs) are obtained for T and cos#b using MC 
simulations for signal and background. The PDFs are then combined to form 
a likelihood ratio C = C S /(C S + C q g), where C s and C q g are the product of 
the signal and background PDFs, respectively. Signal events are selected by 
requiring that £ > 0.9, which has a signal efficiency of about 40% for modes 
with a 7t° and 35% for the all charged particle mode. The systematic error in 
the likelihood ratio requirement is determined using a sample of B + — > D 7r + 
decays. By comparing the yields in data and MC after the likelihood ratio 
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requirement, the systematic errors are determined to be 4% for the modes 
with a 7T° and 6% for the mode without. 

The final variable used for continuum suppression is the p helicity angle, Oh, 
defined as the angle between the direction of the decay pion from the p in the 
p rest frame and the p in the B rest frame. The requirement of | cos 0^1 > 0.3 
is made independently of the likelihood ratio as it is effective in suppressing 
the background from B decays as well as the qq continuum. 

Although continuum events are the largest background, the mode with all 
charged particles has other backgrounds in the signal region from B meson 
decays. To study this, a large MC sample of BB events is used [14]. The 
largest component of this background is found to come from decays of the 
type B — > Dn; when the D meson decays via D — > 7r + 7r~, events can directly 
reach the signal region while the decay D — > K~n + can reach the signal 
region with the kaon misidentified as a pion. Decays with J/ip and ip(2S) 
mesons can also populate the signal region if both the daughter leptons are 
misidentified as pions. These events are excluded by making requirements 
on the invariant mass of the intermediate particles: |M(7r + 7r~) — Mjjo\ > 0.14 
GeV/c 2 , \M(tt + tt°)-M d+ \ > 0.05 GeV/c 2 , \M(it + it-) -M J/t \ > 0.07 GeV/c 2 
and \M(n + 7Y~) — M^2S)\ > 0.05 GeV/c 2 . The widest cut is made around the 
D° mass to account for the mass shift due to misidentifying the kaons in D° 
decays as pions. 

Fig. 1 shows the AE and M bc distributions for the three modes analysed after 
all the selection criteria have been applied. The AE and M bc plots are shown 
for events that lie within 3a of the nominal M bc and AE values, respectively. 
The signal yields are obtained by performing maximum likelihood fits, each 
using a single signal function and one or more background functions. The 
signal functions are obtained from the MC and adjusted based on comparisons 
of B + — > D°7T + decays in the data and MC. All parameters in the signal 
function are fixed except the overall normalization, which is allowed to float. 
The final branching fractions are extracted using the AE signal yield, as this 
provides a kinematic separation of the background modes. In calculating the 
branching fractions, the decay rates of the T(AS) into B + B~ and B°B° are 
assumed to be equal. 

The M bc distribution for all modes is fitted with a single Gaussian and an 
ARGUS background function [15]. The normalization of the ARGUS function 
is left to float and shape of the function is fixed from the AE sideband: —0.25 
GeV < AE < -0.08 GeV and 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 . For the 
mode with only charged pions in the final state, the AE distribution is fitted 
with a single Gaussian for the signal and a linear function with fixed shape for 
the continuum background. The normalization of the linear function is left to 
float and the slope is fixed from the M bc sideband, 5.2 GeV/c 2 < M bc < 5.26 
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GeV/c 2 , \AE\ < 0.3 GeV. There are also other rare B decays that are expected 
to contaminate the AE distribution. For the mode without a 7r°, these modes 
are of the type B° — > h + h~ (where h denotes a n or X), _B — > pp (including 
all combinations of charged and neutral p mesons, where the polarizations of 
the p mesons are assumed to be longitudinal) and B — > Kim (including the 
decays B + -> p°K + , B + -> K* 7r+, 5+ -> K*(1430)°7r + , 5+ -> f (980)K+ 
and _B + — > /o(1370).fC~) [16]. These background modes are accounted for by 
using smoothed histograms whose shapes have been determined by combining 
MC distributions. The three B — > pp modes are combined into one histogram. 
The normalization of this component is allowed to float in the fit due to the 
uncertainty in the branching fractions of the B — > pp modes. Likewise, the 
B — > hh and all the B — > Knir modes are combined to form one hh and one 
Kim component. The normalizations of these components are fixed to their 
expected yields, which are calculated using efficiencies determined by MC and 
branching fractions measured by previous Belle analyses [16,17]. 

The AE fits for the modes with a ir° in the final state have the signal compo- 
nent modelled by a Crystal Ball function [18] to account for the asymmetry in 
the AE distribution. As for the B + — > p°ir + mode, the continuum background 
is modelled by a linear function with fixed slope. Unlike the B + — > p°ir + mode, 
a component is included for the background from the b — > c transition. The pa- 
rameterization for rare B decays includes one component for the B — > Kirn 
modes (B° — > p + K~ and B° — > K* + n~) [19] and one for all the B — > pp 
modes. The normalization of the S — > pp component is left to float while the 
other components from B decays are fixed to their expected yields. 

Table 1 summarizes the results of the AE fits, showing the number of events, 
signal yields, reconstruction efficiencies, statistical significance and branching 
fractions or upper limits for each fit. The statistical significance is defined as 
\J —2 \n(£o/£ max ), where C max denotes the likelihood at the nominal signal 
yield and C is the likelihood with the signal yield fixed to zero. For the p°ir + 
and p + 7r~ modes, the yields obtained from the AE fit are 24.3±@;2 and 44.6li§;|, 
respectively. Both signals have a significance of over 3a and consistent signal 
yields are obtained from the M bc distributions. No excess above background 
is seen in the AE fit for the p°7r° mode and an excess of 14 ± 6 events is seen 
in the M bc distribution. The M bc yield has a significance of 2.3<r and may be 
due to background from rare B decays. An upper limit (UL) is obtained for 
this mode based on the AE yield. The upper limit is calculated using the 
90% confidence level UL of the signal yield (N s ), obtained by integrating the 
likelihood function to be f£ L C(N s )dN s / / °° C(N s )dN s = 0.9, where C(N S ) 
denotes the maximum likelihood with the signal yield fixed at N s . 

The systematic errors in the branching fractions are obtained by quadratically 
summing the systematic uncertainties in the signal yield, tracking, particle 
identification, n° reconstruction and the likelihood ratio requirement. The 
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systematic error in the fitted signal yield is estimated by independently varying 
each fixed parameter in the fit by la. The final results are B(B + — > p°ir + ) = 
(8-Oll:g±8:;5 x 1(T 6 and B(B° -> P V) = (20.8±g:g±i:f ) x 10- 6 where the 
first error is statistical and the second is systematic. For the p°7r° mode, one 
standard deviation of the systematic error is added to the statistical limit to 
obtain a conservative upper limit at 90% confidence of 5.3 x 10~ 6 . 

The possibility of a nonresonant B — > unit background is also examined. To 
check for this type of background, the Mf, c and AE yields are determined for 
different nn invariant mass bins. By fitting the Mb c distribution in tttc invariant 
mass bins with B — > pn and B — > mm MC distributions, the nonresonant 
contribution is found to be below 4%. To account for this possible background, 
errors 3.7% and 3.2% are added in quadrature to the systematic errors of the 
p + 7r~ and p°iT + modes, respectively. The tttt invariant mass distributions are 
shown in Fig. 2. Two plots are shown for the p + n~ and p°7v + modes, one 
with events from the M^ c sideband superimposed over the events from the 
signal region (upper) and one with events from signal MC superimposed over 
events from the signal region with the sideband subtracted (lower). Fig. 3 
shows the distribution of the helicity variable, cosO^, for the two modes with 
all selection criteria applied except the helicity condition. Events from pir 
decays are expected to follow a cos 2 9 distribution while nonresonant and other 
background decays have an approximately uniform distribution. The helicity 
plots are obtained by fitting the M^ c distribution in eight helicity bins ranging 
from —1 to 1. The M^ c yield is then plotted against the helicity bin for each 
mode and the expected MC signal distributions are superimposed. Both the 
7T7T mass spectrum and the helicity distributions provide evidence that the 
signal events are consistent with being from p7r decays. 

The results obtained here can be used to calculate the ratio of branching frac- 
tions R = B(B° -> p ± ir^)/B(B + -> p°7r+), which gives R = 2.6 ± 1.0 ± 0.4, 
where the first error is statistical and second is systematic. This is consistent 
with values obtained by CLEO [20] and BaBar [21,22] as shown in Table 2. 
Theoretical calculations done at tree level assuming the factorization approx- 
imation for the hadronic matrix elements give R ~ 6 [3]. Calculations that 
include penguin contributions, off-shell B* excited states or additional im res- 
onances [4-8] might yield better agreement with the the measured value of 
R. 

In conclusion, statistically significant signals have been observed in the B — > 
p-n modes using a 31.9 x 10 6 BB event data sample collected with the Belle de- 
tector. The branching fractions obtained are B{B + — > p Q n + ) = (S.Ol^o-a?) x 
10" 6 and B(B° -> p ± vr T ) = (20.81^11?) x 10~ 6 where the first error is statisti- 
cal and the second is systematic. No evidence was seen for the mode B° — > p°7r° 
and an upper limit of 5.3 x 10~ 6 at 90% confidence level was obtained. 
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Table 1 

AE fit results. Shown for each mode are the number of events in the fit, the signal 
yield, the reconstruction efficiency, the branching fraction (B) or 90% confidence 
level upper limit (UL) and the statistical significance of the fit. The first error in 
the branching fraction is statistical, the second is systematic. 



Channel 


Events 


Signal Yield 


Eff.(%) 


B/UL( xlO" 6 ) 


Significance 


A+ 
A° 


154 
301 
116 


24 3+ 6 9 

44 R+ 12 - 8 

-4.4 ±8.5 


9.6 
6.8 
8.5 


o n+2.3+0.7 
°' u -2. 0-0.7 

90 Q+6. 0+2.8 
zu -°-6. 3-3.1 

< 5.3 


4.4,7 
3.7(7 



Table 2 

A comparison between the results of the Belle, BaBar and CLEO experiments for 
the B — > pir branching fractions and the ratio R = B(B° — > p ± ir Zf )/B(B Jr — > p°7r + ). 
The CLEO results come from reference [20]. The BaBar results for p + i:~ and 
come from references [21] and [22], respectively. 



Experiment 


B(B° p^T) x 10- 6 


B{B+ -> p°vr+) x 10" 6 


R 


Belle 


on c+6.0+2.8 
zu -°-6.3-3.1 


o n+2.3+0.7 
°- u -2.0-0.7 


2.6 ± 1.1 


BaBar 


28.9 ±5.4 ±4.3 


24 ± 8 ± 3 


1.2 ±0.5 


CLEO 


27.6±|;| ± 4.2 


10.41^ ± 2.1 


2.7 ±1.3 
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Fig. 1. The AE (left) and M^ c (right) fits to the three B — > pn modes: p°7r + , 
p + ir~ and p 7r°. The histograms show the data, the solid lines show the total fit 
and the dashed lines show the continuum component. In (a) the contribution from 
the B — > pp and B — > hh modes is shown by the cross hatched component. In (c) 
the cross hatched component shows the contribution from the b — > c transition and 
B —> pp modes. 
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Fig. 2. The M(inr) distributions for B° -> p ± vr T (left) and 5+ 
events in the signal region. Plots (a) and (b) show sideband events superimposed; 
plots (c) and (d) show the sideband subtracted plots with signal MC superimposed. 
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Fig. 3. The p meson helicity distributions for B° — > p^i^ (a) and -B + 
Signal MC distributions are shown superimposed. 
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